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Abstract 

Superpositions of different statistics on different time or spatial scales (in short, 
superstatistics) can naturally lead to an effective description by nonextensive sta- 
tistical mechanics. We first discuss the role of escort distributions within the super- 
statistical framework, and then briefly describe recent physical applications of this 
concept to turbulent and patttern forming systems. 
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1 Introduction 



In this paper we will briefly review the so-called superstatistics concept [1]. 
It yields a simple and plausible argument why nonextensive methods, as de- 
veloped by Tsallis and many others [2,3,4], are relevant for nonequilibrium 
systems with large-scale spatio-temporal fluctuations of an intensive parame- 
ter. 

Superstatistics is a superposition of two (or even more) different statistics: One 
given by ordinary Boltzmann factors, and another one given by large-scale fluc- 
tuations of one (or several) intensive parameters (e.g. the inverse temperature). 
The corresponding stationary probability distributions arise as a convolution 
of the various statistics. Tsallis statistics is a special case that arises naturally 
in this approach for a x 2 -distributed parameter, if the marginal distributions 
are formed. But other generalized statistics are possible as well. Recently it has 
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become clear that the concept is not only a theoretical construct but of prac- 
tical physical relevance. For example, it is possible to develop superstatistical 
models that well describe the statistical properties of measured accelerations 
in Lagrangian turbulence experiments [5,6], or the statistics of defect velocities 
in pattern forming systems [7] , or the measured energy distribution of cosmic 
rays [8,9]. 



2 Basic idea of superstatistics 



Consider a driven nonequilibrium systems with spatio-temporal fluctuations of 
an intensive parameter (3. These fluctuations are externally produced, by con- 
stantly putting energy into the system which is dissipated. The intensive pa- 
rameter may be the inverse temperature, or an effective chemical potential, or 
a function of the fluctuating energy dissipation in the flow (for the turbulence 
application). Locally, i.e. in spatial regions (cells) where (3 is approximately 
constant, the system may be described by ordinary statistical mechanics, i.e. 
ordinary Boltzmann factors e _/3£; , where E is an effective energy in each cell. 
In the long-term run, the system is described by a spatio-temporal average 
over the fluctuating (3. 

Suppose the probability density of (3 is f{(3). By integrating over (3 one obtains 
marginal distributions of the form 

oo 

m = J f(P)^f )P (E)e^ E df3, (1) 



where p(E) is the density of states and Z(f3) is the normalization constant 
of p(E)e~ l3E for a given (3. Eq. (1) describes the long-term behaviour of the 
system (we have chosen type-B superstatistics [1]). 

Let us consider three important examples (many more are possible). 

1. (3 is distributed according to the ^-distribution 




n is a parameter. The average of the fluctuating (3 is given by 

oo 

{(3) = J (3f((3)d(3 = fa (3) 



2 



and the second moment of (3 is given by 

f-4 



Let E = |m 2 be a kinetic energy. The marginal distribution (1), obtained by 
integration over all /3, yields the generalized canonical distributions of nonex- 
tensive statistical mechanics [2,3,4] 

1 1 

p(«) = 7- rzr: (5) 



where Z q is a normalization constant and 



«= 1 + ;T7T < 6 > 

(7) 

The distributions (5) maximize the Tsallis entropies S g subject to suitable 
constraints [2,3,4]. 



2. p is log-normally distributed: 
1 [ -(log 



m 1 expl^} (8 ) 



The average /?o of the above log-normal distribution is given by /3q = m^/w 
and the variance by a 2 = m 2 w(w — 1), where w := e s . One obtains the 
superstatistics distribution 

p(u) = JL /T 1/2 exp{ ~ (1 ^ )2 | e-^« 2 . (9) 
o ^ ^ 



The integral cannot be evaluated in closed form, but it can be easily numeri- 
cally evaluated and compared with experimental data [5,6,10]. 

3. Yet another interesting example of a superstatistics is obtained if we assume 
that the temperature T = (3~ l (rather than (3 itself) is x 2 -distributed. In this 
case the corresponding marginal distribution p(u) has exponential tails [11]. 
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3 Escort distributions in superstatistical systems 



For a given probability distribution {p^ of microstates i and a given parameter 
q the escort distribution {Pi} is denned as [12] 



Both types of distributions, {pi} and {Pi}, notoriously occur in the formulation 
of nonextensive statistical mechanics [3,4,13], like two spin-degrees of freedom 
which are not there in classical statistical mechanics (q = 1). We will now 
show that within the superstatistical approach these two degrees of freedom 
can be easily understood. 

Generally, it makes sense to /^-average typical thermodynamic relations if the 
local cells are sufficiently large and if they are in local equilibrium. In this case 
one has for each local cell 

S = j3(U-F), (11) 



where S is the ordinary Boltzmann-Gibbs-Shannon entropy, U is the internal 
energy and F the free energy associated with each cell. The averaged version 
reads 

oo oo oo 

J f((3)Sd(3 = J (3f{(3)Ud(3 - J j3f(p)Fdp. (12) 





Note that on the right-hand side we obtain averages formed with /3f(/3). For 
X 2 -distributed (3 this is like transforming n — > n + 2 in the distribution (2). 
Since 

^T = ^ (13) 



and 



n + 3 



(14) 



this means one is forming averages with respect to the escort distributions on 
the right-hand side of eq. (12), as required in the modern versions of nonex- 
tensive statistical mechanics [3]. 
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Our argument is valid for arbitrary energies E. We see that within the super- 
statistical /^-averaging approach entropic quantities most naturally correspond 
to the original distribution p iy energetic quantities such as U or F most natu- 
rally to the escort distribution P, h since there is an additional factor (3 on the 
right-hand side of eq. (12). 

In general, the averaged Shannon entropy on the left-hand side does not co- 
incide with the Tsallis entropy S q , which satisfies [2,3,4] 

S q = 0(U q -F q ). (15) 



However, it might be interesting to construct special situations where eq. (12) 
and (15) are equivalent. 



4 Applications 



Langevin equations where the parameters fluctuate on a large spatio-temporal 
scale yield a concrete dynamical realization of superstatistics [14]. For example, 
to model the motion of a single test particle in a turbulent flow one considers 
a superstatistical extension of the Sawford model [15]. The Sawford model 
describes the joint stochastic process (a(t),u(t),x(t)) of acceleration, velocity 
and position of a Lagrangian test particle in a turbulent flow by the stochastic 
differential equation 



+ ^ 2 u(T L 1 + t^)T L H^L(t) (16) 
ii = a (17) 
x = u, (18) 

where L(t) is Gaussian white noise. Note that in this model the acceleration a 
and not the velocity u is driven by white noise, so the meaning of the variables 
is different as compared to an ordinary Brownian particle. Tl and t rj are two 
time scales, with T L » t v and 



2a v l l 2 
Coe 1 / 2 ' 



(19) 



(20) 
(21) 
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Fig. 1 Histogram of accelerations a as measured in [16] and the log-normal 
superstatistics prediction eq. (9) with s 2 = 3.0. For comparison, the figure also 
shows Tsallis statistics (eq. (5) with q — 1.5). 

e is the energy dissipation, v the kinematic viscosity, Co,ao are Lagrangian 
structure function constants, and o\ is the variance of the velocity distribution. 
The Sawford model with constant coefficients predicts Gaussian stationary 
distributions for o and u, and is thus at variance with recent measurements 
[16]. However, a superstatistical generalization of the Sawford model [5,6] well 
fits the data. 

For Tl — > oo one can derive [6] that the fluctuating parameter (3 is given by 



i.e. the fluctuating energy dissipation e is proportional to (3~ 2 ^. Fig. 1 shows 
the measured probability density of the acceleration of a Lagrangian test par- 
ticle in a turbulent flow as obtained in the experiment of Bodenschatz et al. 
[16]. Log-normal superstatistics with s 2 = 3.0 yields an excellent fit [6]. 

Similar superstatistical models can be formulated for pattern forming systems, 
just that the meaning of the variables in the Langevin equation is different. The 
probability density of defect velocities in inclined layer convection experiments 
can be quite precisely measured. As shown in Fig. 2, it quite precisely coincides 
with a Tsallis distribution with q ps 1.46 [7]. A superstatistical model well 
describes various stochastic properties of these defects. 

Superstatistical models also have applications in high energy physics, for ex- 
ample as simple models to explain the measured statistics of cosmic rays [8,9]. 




(22) 
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Fig. 2 Measured distribution of positive and negative defect velocities in in- 
clined layer convection and comparison with a nonextensive canonical distri- 
bution of type (5) with q — 1.46 (more details in [7]). 
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